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SUMMARY 


A study to determine the effects of various error sources on a decoupled 
control system was performed in connection with the longitudinal control system 
on a simulated externally blown jet-flap STOL aircraft. The system employed 
the throttle, horizontal tail, and flaps to decouple the forward velocity, 
pitch angle, and flight-path angle. The errors considered were (1) imperfect 
knowledge of airplane aerodynamic and control characteristics, (2) imperfect 
measurements of airplane state variables, (3) change in flight conditions, and 
(4) lag in the airplane controls and in the engine response. The error due to 
change in flight conditions results from the use of constant gains to represent 
a range of flight conditions. 

Despite the large errors employed in the study, the effects of the various 
errors on the decoupling process were generally minor. Because changes in 
flight conditions had little effect, constant gains in a decoupled control sys- 
tem for certain flight regimes would be feasible, thereby avoiding the need for 
an onboard computer. 

The largest effects on decoupling were caused by control lag and by change 
in flight condition during speed-command maneuvers wherein significant coupling 
in the flight-path angle was produced. However, these effects could be elim- 
inated by retrimming the flight-path angle when large speed changes were 
required, such as in high decelerating approaches. In addition, lag effects 
could be eliminated by including the lag in the design of the decoupled control 
system. 

Except for the effects on the flight-path angle during speed-command maneu- 
vers, the largest effect of any error source generally occurred in the commanded 
response quantity; the other responses were materially unaffected, as were the 
control requirements. It was shown that if a pilot could recognize these dis- 
crepancies in the command response, he could compensate with an additional com- 
mand input, with only small increased effects on the decoupling process. 


INTRODUCTION 

Decoupling is one of the more sophisticated forms of automatic control 
devised in recent years. (See, for example, refs. 1 to 4.) The method of air- 
plane decoupled control can be briefly stated in three steps: 

(1) Choose the airplane response quantities to be controlled. 

(2) Provide the pilot with one input channel for each quantity. 

(3) Design the control system so that an input command in one channel 
gives desirable response to the command input and no response in the other 
channels . 



An automatic control system of this type would have many practical appli- 
) cations as, for example, use in strafing fast-moving targets where constant 
pitch angle must be maintained during changes in flight-path angle and speed. 
Decoupled control would also be highly desirable for STOL aircraft operation 
where severe control requirements may exist during approach and landing because 
of excessive coupling between the longitudinal flight variables, flight-path 
angle, forward speed, and pitch attitude. The control problem is further com- 
plicated for this powered lift airplane because a sizable flight-path-angle 
response is associated with the throttle; furthermore, if jet engines are the 
power plants, a considerable time lag occurs between the throttle movement and 
the change in engine thrust. 

Application of decoupled control systems to certain types of aircraft has 
been studied extensively over the past few years. (See, for example, refs. 5 
to 8.) Knowledge of the effects of various error sources is necessary for 
future practical applications. Therefore, this paper investigates the magni- 
tude of these effects on the decoupling process. The effects studied include 
the effect of errors in the (1) knowledge of the airplane aerodynamic and con- 
trol characteristics, (2) measurements which operate on the feedback gains, 

(3) use of constant gains when the flight conditions are actually changing, 
and (4) lag in the control actuators. 

The present study was conducted using a simulated STOL aircraft for decou- 
pled longitudinal control in a landing approach condition. The computer pro- 
gram reported in reference 9 was used for the decoupling calculations. This 
program provides complete decoupling because both the transient and steady-state 
airplane responses are decoupled. In addition to the decoupling, the program 
permits the choice of the zeros and poles of each input-output transfer function 
by specifying the polynomial coefficients of the numerator and denominator. 


SYMBOLS 


A matrix defined by aircraft stability coefficients 

/V 

A matrix = A + BF 
a element of A matrix 

B matrix defined by aircraft-control coefficients 

b element of B matrix 

C matrix relating desired output vector to state vector 

Cj[, pitching-moment coefficient 
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Cx 


aircraft weight in coefficient form 
longitudinal-force coefficient 
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mean aerodynamic chord, m 

matrix of feedback gains used in decoupled controller (see appendix A) 
element of F matrix 

matrix of feedforward gains used in decoupled controller (see 
appendix A) 

acceleration due to gravity, m/sec^ 
identity matrix 

moment of inertia about Y-axis, kg-m^ 
indices 

mass of airplane, kg 
pitch rate, rad/sec 
wing area, m^ 

Laplace operator 
time, sec 

velocity component along X-axis, m/sec 
m-vector of control variables 
equilibrium or trim airspeed, m/sec 
m-vector of commanded pilot inputs 
body-axis system 
n-vector of state variables 

m-vector of state variables to be controlled in decoupled manner 
angle of attack, deg 
flight-path angle, deg 

incremental change from equilibrium or trim condition 
deflection of trailing-edge flaps, deg 
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horizontal-tail deflection, deg 
throttle deflection 

short-period damping ratio for closed-loop (decoupled) airplane 

damping ratio for horizontal- tail servo 

damping ratio in speed response 

damping ratio in y response 

pitch angle, deg 

air density, kg/m^ 

flap response time constant, sec 

horizontal-tail response time constant, sec 

thrust response time constant, sec 

speed response time constant, sec 

flight-path-angle response time constant, sec 

longitudinal short-period natural (undamped) frequency for closed-loop 
(decoupled) airplane, rad/sec 

natural frequency of horizontal-tail servo, rad/sec 

natural frequency in speed response, rad/sec 

natural frequency in y response, rad/sec 


Aircraft stability and control coefficients: 


3Cx 9Cz 9Cjn 

^^•Sf ’ 36f ~ 36f ^<Sf ~ 36f 

3Cx 3Cz 3Cm 

^^•Sth " 36th " 36th ” 36th 

3Cx 3Cz 3Cjn 

^^<St 95^ ^^6t " 96^. *^6t ■ 96t 
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Superscripts; 

-1 matrix inverse 

' perturbations from equilibrium 

Subscr ipts ; 

0 equilibrium or trim condition 

c commanded by pilot 

Dot over symbol denotes differentiation with respect to time. 


DESCRIPTION OF SIMULATED AIRPLANE 

The STOL aircraft simulated in this study is the clustered-engine airplane 
aerodynamically described in references 10, 11, and 12. The airplane (fig. 1) 
is a high-wing jet transport with four high bypass ratio turbofan engines. The 
four engines yield a maximum total thrust of 147 057 N. The pertinent physical 
characteristics of the simulated aircraft including maximum control-surface 
deflection and deflection rate are presented in table I. The aerodynamic char- 
acteristics are presented in table II. Values are shown for two flight condi- 
tions, a = 10° and a = 5°. In both cases, the nominal speed is 30.48 m/sec 
and the thrust coefficient is 1.87 (one-half power). Where a = 10°, two values 
are changed which account for thrust effects on Cg and that were not 

included in the standard set of coefficients. 
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DECOUPLED CONTROL 


The approach taken in providing independent longitudinal control of pitch 
angle, flight-path angle, and forward velocity is depicted in sketch (a) . 
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equations of motion 
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Sketch (a) 


The feedforward gain matrix G operates on the pilot command inputs in 
order to provide output solely in the commanded airplane responses. Both gain 
matrices P and G operate to decouple and yield desired responses. The state, 
variables must be continuously measured and operated on by the feedback gain 
matrix F to move the airplane controls automatically in order to maintain the 
desired decoupled responses. The C matrix allows any combination of the state 
variables to be selected as the output quantities to be decoupled. Figure 2 
illustrates the effectiveness of decoupled control. The results in figure 2(a) 
pertain to a horizontal-tail step input of about 2° for the unaugmented (no 
decoupled control) system. Not only is the coupling between the three response 
quantities apparent, but the intended response in 6 is sluggish. For the 
decoupled system, figure 2(b) shows how the controls automatically move to main- 
tain the u and y responses at 0 for a pitch command of about 3°. Also, the 
0 response has been improved by selecting (Oj^ = 2 rad/sec and C = 0.7 as the 
closed-loop (decoupled) dynamic characteristics. 

The most sophisticated implementation of decoupled control would employ an 
onboard computer to determine the time-varying adaptive gain matrices F and G. 
However, this degree of complexity may not be required if the decoupled control 
is limited to certain flight regimes, such as the approach and landing phase. 

In the present study these F and G gain matrices are treated as constants, 
thus eliminating the need for an onboard computer. There are various methods of 
determining the required F and G gains for longitudinal decoupled control. 

One method, which both simulation and flight testing (ref. 8) have proven to be 
practical, is devised in reference 7 for decoupled control of the steady-state 
conditions. Other methods have been developed for complete decoupling in order 
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to include the transient conditions. The method of determining F and G in 
reference 3, along with the computer program in reference 9, is used for the 
present study. 

A general development of the method for complete decoupling is given in 
appendix A. Because of its clarity and simplicity, Rekasius' method (ref. 1) 
for decoupling analysis was chosen instead of the abstract geometric and alge- 
braic methods of most decoupling analyses. 


METHOD OF ANALYSIS 

All cases in the present analysis pertain to decoupled control of the lon- 
gitudinal response quantities u, 6, and Y* The values of the coefficients 
for the A and B matrices, as determined from equations (Al) to (A3) and 
from table II, are presented in table III. Unless otherwise noted, the results 
in this paper pertain to a = 10° (see table 111(a)) which has been denoted 
the standard case for this analysis. 

The analysis presented in this paper is quantitative in nature and is 
intended to give an overall indication of the effects of various error sources 
on the decoupling process. The errors considered were (1) imperfect knowledge 
of airplane dynamics, (2) imperfect measurements of airplane state variables, 

(3) change in flight conditions, and (4) lag in the airplane controls and in 
the engine response. 

The results of the study were obtained from computer time-history simula- 
tions of step-command inputs using values of the F and G gain matrices 
obtained from the computer program of reference 9. Table IV presents typical 
values of the gain matrices and the airplane design response characteristics 
used in the analysis. 


DISCUSSION OF RESULTS 

The results of the investigation are presented primarily as computer- 
simulated time histories of airplane responses and control motions. Time his- 
tories pertaining to various types of errors are compared with the nominal (no 
error) time histories to show the magnitude of the effects of the errors 
imposed on the system. 


Effect of Imperfect Knowledge of Airplane Characteristics 

In order to compute the constant gains for the decoupled control system 
designers need the airplane aerodynamic and control characteristics, as 
reflected by the A and B matrices (appendix A) . These characteristics are 
determined from the airplane stability and control derivatives which must be 
obtained from wind-tunnel tests, flight tests, or theoretical calculations. 
There is always some question about the accuracy of these values and, hence, 
about the computed gains required for decoupling. Nonlinear effects also 
account for error in the A and B matrices. This error is introduced by 
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linearizing the airplane equations of motion in the decoupling calculations 
(appendix A) . 

Effects on airplane decoupled responses .- Tables V and VI are presented 
as a quantitative summary of the effects of errors in the different A and 
B elements on decoupling. In each table the stability or control derivatives 
mainly responsible for each element are given in the different blocks of the A 
or B matrix found at the top of the tables. As noted by the zero b 2^1 and 
^3,1 elements in table VI, the results apply to the standard case. 

For the A matrix, the elements were reduced to zero, one at a time, to 
determine their individual effects on decoupling. Decoupling was more sensitive 
to errors in the B-matrix elements; hence, the elements were reduced to one-half 
their nominal values only. In all cases the values of the F and G matrices 
remained the same. The results of changing the values of the elements are shown 
in the charts at the bottom of the tables. For example, for an error in Si2 2 > 
the chart shows that decoupling is not lost for a pitch command, but that an 
overshoot of 20 percent in the commanded pitch value occurs. (The terms "over- 
shoot" and "undershoot" refer to the ensuing steady-state error during a command 
maneuver.) The results for pitch command are similar for an error in 32 3* 

This overshoot effect is not serious, as the pilot can compensate for this 
effect, as shown in subsequent examples. For an error in 33 3 , the chart indi- 
cates that decoupling is lost in the flight-path-angle response for a pitch, com- 
mand; decoupling is not significantly affected for the other two commands. For 
an error in the elements in the fourth column of the A matrix, the charts show 
no significant effect on decoupling for either a 6 or y command, and only 
small effects for a speed command. For the blocks showing partial decoupling, 
the effects are not considered serious because of the severe errors that were 
imposed on the A and B elements. 




Further tests vrfiose results are not presented in the tables were made with 
^ and Cni^ ^ 0. Changing b 3^1 ^€ 2 ^ to one-half its nominal value had 


no significant effect on decoupling. A similar change in b 2 ^ 1 , however, 
destroyed decoupling in y for 6 and u commands and caused a large 
(70-percent) overshoot in y for a y command. Effects on the control 
motions were not excessive in any case. 


Effects on control motions .- Some effects of the A- and B-matrix errors 
on the control motions, as well as on the Commanded responses, are shown in 
figures 3 to 5. The effects on the controls may be important because of con- 
trol displacement or rate limits. In regard to A-matrix errors (fig. 3), the 
flap and throttle are the controls most affected. The maximum flap deviation 
from the nominal is about 15 percent of full deflection, and the maximum devia- 
tion from the nominal for the throttle is about 20 percent. (Fig. 3(b) shows 
a small overshoot in y for an error in a 2 ^ 3 , not noted in table V because 
decoupling in 0 has been completely lost.) B-matrix errors which affect 
decoupling show little effect on the controls for a 6 command (fig. 4(a)). 

For a y command, an error in b 4 ^ 3 ^ is apparently the most serious with 
regard to the controls; the steady-state throttle deflection increases over 
150 percent of nominal and would exceed its limit value. Figure 5(a) is 
presented to show the amount of coupling in 0 for an error in b 2 ^ 2 - 


8 



commanded value of Yf however, is still achieved (Y =9 - ot) . For complete 
decoupling, y = • 

As shown in figure 5(a), an error of 1° is obtained in 6 . If the pilot 
could recognize this error, he could compensate for it by employing a 9 com- 
mand of 1.12° as shown in figure 5(b). The additional 12-percent command is 
needed to account for the undershoot due to the b2^2 error. (See table VI.) 
Also, for a pitch command, table VI shows that the b2^2 error produces no 
significant effect on Y • By using the 9 command in figure 5(b), the 9 
error is reduced to zero, with the command value of 6 ° still being maintained 
for Y • 

Summary of effects .- The overall effects of errors in the A and B 
matrices can be summarized as follows: The more important A-matrix elements, 
as far as sensitivity to decoupling is concerned, appear to be 84,1 (C^) and 
those in the third column (a derivatives). For the B matrix, errors in 
^ 4,1 ^ *^4,3 have a coupling effect on u for 9 and Y 

commands, while the b2^2 element essentially affects only 9 during a 

Y command. Also, for the case where the nominal value of b 2 ,l was not zero, 
errors in this element produced coupling in Y for 9 and u commands. 

For A-matrix errors, effects on decoupling as well as changes in the con- 
trol motions were essentially linearly related to the magnitude of the errors. 
This linearity was not the case for B-matrix errors because of their highly 
nonlinear effect on the decoupling process. In some instances, changing the B 
elements to zero magnitude prevented decoupling , whereas the one-half nominal 
magnitude error had only a marginal effect. These results indicate that a 
better knowledge of the B matrix than of the A matrix is required in 
designing adequate decoupled control systems. 

Effect of Imperfect Measurements of State Variables 

In the decoupling process the feedback gain matrix F operates on the 
measured values of the airplane state variables 9, 9, a, and u. In any 
practical' design, errors in these measurements can have an effect on decoupling 
In extreme cases these errors can become quite large owing to instrument mal- 
function. The measurement-error effects were investigated, both individually 
and in combination. For the present analysis, these measurement errors were 
simulated by changing the values in certain columns of the F and BF matri- 
ces by various amounts. Equations (A12) and (A14) are repeated here to show 
that these changes are equivalent to prescribing errors in the various measured 
response quantities, x. 

u = Fx + Gv 

X = (A + BF) X + BGv 

For example, according to the matrix which indicates the arrangement of the 
equations of motion as given in table III, a 20-percent increase in the values 
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of the third columns of F and BF would represent a 20-percent error in the 
a measurement. 

Individual effects on decoupling .- Examples of the effects of errors in 
the various measurements are presented as time-history comparisons in figures 6 
to 9. These effects are summarized in table VIJ. The table clearly indicates 
that measurement inaccuracy in the pitch rate 0 is not important. Also, the 
table shows that the dominant effect of the errors in each of the other three 
variables occurs in the particular commanded response quantity in which the 
error exists. For example, with measurement error in a (which corresponds to 
-y in the decoupled system) , a large effect on y and small effects on 0 
and u are noted for a y-command maneuver and little or no effect on 0, y, 
and u for either a 0 or u command. Time-history examples of these effects 
are presented in figure 6. For the cases in figure 6 the effects in the air- 
plane responses and control motions appear to be essentially linear with respect 
to the magnitude of a-measurement error. 

With regard to velocity-feedback error, table VII shows that this error is 
important only in a u-command maneuver. An example time history of this type of 
maneuver is shown in figure 7. The figure indicates that the control motions 
are not significantly affected. 

The time histories in figure 8 illustrate the effect of 0-measurement error 
for a 0-command maneuver. The major effect of the error is in the 0 response. 
Two nominal cases are included in the figure to show how the controls react for 
two different sets of decoupled responses. In the one case u, 0, and y are 
decoupled (controlling y) , and in the other, u, 0, and a are decoupled 
(controlling a) . In order to avoid confusion, nominal values of y and a 
are not plotted in figure 8. For the controlled y case, nominal a, of 
course, equals 0 and in the controlled a case, nominal y also equals 0. 

It is shown that the measurement error does not affect the a response; how- 
ever, there is a small effect in y inasmuch as y depends on 0 in the 
equations of motion. 

Combined effects on decoupling .- The effects of combined measurement errors 
are illustrated in table VIII. The results in the table show primarily that 
combined errors have about the same effect as individual errors. For example, 
in table VIII (b) the effect of combined errors is the same as the effect of an 
a error only; that is, the effect of adding other errors is insignificant. 

In table VIII (a), for a 0 command, only the 0-measurement error is important. 
The effects of combined measurement errors are not shown for u-command maneuvers 
because they are negligible. Thus, for a u command, only the u-measurement 
error is important. 

Linearity characteristics of error effects .- During the investigation it 
was noticed that the effects of 0-measurement errors were nonlinear. An example 
of this nonlinearity is apparent in table VIII (a) for positive and negative 0 
error. Another example is shown in time-history form in figure 9. It also 
became apparent that the error in the u or y commanded response quantity 
depends on the magnitude of the nominal value of the corresponding element in 
the A matrix. For example, for speed response the pertinent element would be 
®4,4 (according to the arrangement shown in table III), since this element 
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corresponds to the term in the velocity equation (u row) that would be 
affected by u-measurement error. For Y response the pertinent element would 
be 23 , 3 . The linearity characteristics of the error effects were analyzed by 
developing the following equations. 

From the decoupled equations ( (A13) and (A14) ) 

y = Cx 

X = Ax + BFx + BGv 

let yj[ = X|j where X|j is the commanded response. Then the kth augmented 
equation of motion can be written as 


n m n m 

^ akjXj + ^ bkrfrjXj + 21 d) 

j=l r=l j=l r=l 

1=1 


The steady-state error in the response xj^ caused by an error in the mea- 
surement of X|j can be determined as in the following discussion. For steady 
state, xjj = 0. Assuming a relative error E in the measured X|^ gives the 
following expression for the measured quantity: 

*k,meas = (1 + E)xjt 

Substituting this expression into the term in equation (1) containing the feed- 
back matrix F yields 


®kk + (1 + 



*^krf rk 


m 


Xk = - 


r=l 

Z=1 


^kr9rI^Z ” ^kj^^j 


(j k) 


( 2 ) 


where the term aj^jXj represents the summation of all response-quantity terms 
other than the commanded response. Let R represent the right-hand side of 
equation (2) . For no error in measuring xj^, the equation for the nominal value 
of X|^ can be written as 



m 


bkrfrk xk 


r=l 


/ 


n 


= R 


(3) 


where x^^ represents the nominal value. 

m 

Replacing a^k with a and ^ l^kr^rk with bf gives the relative 

r=l 

error in the steady-state response caused by measurement error in Xk only 
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(using eq. (3) for nominal and eq. (2) for actual) : 


- Xkj 

^kn 


a + (1 + E)bf a + bf 
R 


a + bf 


^k “ *kn a + bf 


XI, a + bf + Ebf 

•'n 


- 1 


*k ~ *k^ -Ebf 


a + bf + Ebf 


*k “ ^kp -E 


Xk a 

1 + E + — 
bf 


^k “ ^kp -E 


xif E + c 

'^n 


where 


(4) 


a 

c = 1 + — 
bf 

An example of results determined from equation (4) is presented in fig- 
ure 10. It is of interest to note that the curve for c = 1, which represents 
the maximum possible error, applies to G under any condition inasmuch as the 
pertinent a element is always zero. Examination of figure 10 reveals the 
nonlinear characteristics of the response error with regard to both measurement 
error and the nominal value of the pertinent element of the unaugmented A 
matrix. A noticeable difference in these characteristics exists between posi- 
tive and negative measurement error. Negative measurement errors have the 
largest effect, especially in the case of the 6 response, where c = 1. 
Positive measurement error refers to a bias error in the measurement which pro- 
duces a decrease in the absolute magnitude of the response quantity, and con- 
versely for negative measurement error. (For example, see fig. 9.) 

Because the magnitude of the commanded response error depends on the nomi- 
nal value of a, the error in the y or u commanded response could be reduced 
by altering the A matrix with a change in the aircraft aerodynamic design. 

Compensating for measurement error .- The results of the overall error anal- 
ysis have shown that measurement error in a given response quantity affects the 
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conrmanded response of that particular quantity primarily, with only small 
amounts of coupling in the other response quantities. In general, this effect 
on the commanded response would not be critical in a decoupled control system 
because the pilot would normally recognize the effect on the response he is 
commanding and could compensate for this error. One such case is shown in fig- 
ure 11. The short-dashed line in the figure duplicates the one in figure 8 and 
shows a 28-percent error in the pitch response. In anticipating the error and 
by compensating with an additional degree of pitch command, as represented by 
the long-dashed curve, the pilot can obtain the originally desired pitch angle 
of 3°. The additional effect on the controls and on the coupling of the other 
airplane response quantities is minor. 


Effect of Change in Flight Conditions 

In the design of a decoupled control system the F and G gain matrices 
are determined by the values of the A and B matrices which, in turn, depend 
on the airplane stability and control derivatives. (See appendix A.) These 
derivatives vary with the airplane flight conditions: speed, angle of attack, 

and thrust coefficient, primarily. As previously noted, a decoupled control 
system employing a set of constant gains corresponding to a certain set of 
flight conditions would be advantageous in any practical implementation. This 
section investigates the feasibility of using constant gains to approximate a 
decoupled control system operating under varying flight conditions. The 
results give an indication of the magnitude of the flight condition effects on 
the decoupled airplane responses as well as on the required control motions. 

Effect of speed .- The effect of changes in speed on the A and B 
matrices is shown in table IX. These changes in the elements of A and B 
affect the decoupling process inasmuch as the gains are calculated only for 
nominal speed ug. In order to determine the effects of speed, the A and* B 
elements were changed according to the ratios in table IX (Ot was assumed not 
to change; hence, the derivatives did not change) whereas the F and G gain 
matrices that were calculated for the original nominal speed of 30 m/sec were 
left unchanged. Time histories were then simulated for individual step commands 
in the airplane responses 9, Y, and u. An example is shown in figure 12 and 
the results of the analysis are summarized in table X. 

Effects on airplane decoupled responses .- As shown by the data in table X, 
speed has an insignificant effect on the decoupling process for 6 and Y com- 
mands. The commanded values are generally obtained with little coupling in the 
other response quantities. Some noticeable effect, however, is apparent for the 
speed-command maneuver. The major effects are seen to be the 1.4° coupling in 
Y and, in the case where Cgg ^ and h ^ large error in the com- 

manded u response. (It is of interest to note that a 1.5-m/sec speed change 
for the unaugmented airplane produced a Y change of 4.5°. ) This case had two 
additional nonzero elements in the B matrix (table III) which resulted in dif- 
ferent F and G gain matrices (table IV) . The speed effects on the u-command 
maneuver are quite different for this case, the difference being caused by the 
change in the nominal values for the F gain matrix. The fact that speed 
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effect depends on the nominal F matrix signifies that this effect would vary 
for different types of aircraft. 

It is important to state that the pilot could compensate for the large 
error in the u response during a speed-command maneuver in the same manner 
described in the previous section where measurement error resulted in large 
discrepancies in the commanded responses. 

The pilot could compensate for the large coupling in y in a manner sim- 
ilar to that shown in figure 5(b). Table X shows that a y command for cor- 
recting Y would have no coupling effect on u or 6. In addition, simulator 
studies in reference 7 show that excessive coupling in y during large speed 
changes (high decelerating approaches) was eliminated by retrimming the flight- 
path angle, either automatically or by the pilot. 

Effects on control motions .- The results of flight-condition change on the 
control motions are not included in table X. There were only slight changes 
noted in the control motions for all cases except the speed-command maneuver. 

A sample case is shown in figure 12 vrtiere the speed has been changed from 
30.48 m/sec to 60.96 m/sec. The speed has the largest effect on the flap, 
vrfiich is a change of about 30 percent from the nominal value. Of course, in 
all cases the controls as well as the response quantities exhibit a noticeable 
change from nominal in the time response as illustrated in figure 12. This 
change is caused by the fact that the & matrix has been altered by changing 
the A-matrix stability derivatives, which is equivalent to changing the aug- 
mented airplane dynamic characteristics. 

Effect of angle of attack .- The effects of changing the angle of attack a 
from 10° to 5° are shown in figure 13 for step commands in 6, y, and u. 
Changes in a affect the airplane stability and control derivatives vrtiich, in 
turn, affect the A and B matrices. (See tables II and III.) As shown in 
table III, a change of 5° in a produces changes of about 5 to 10 percent in 
the A and B coefficients. The effects of the a change on decoupling were 
investigated in the same manner as the speed effects. 

The time-history examples in figure 13 show that an a change had only 
minor effects on the decoupling process. In all three cases the correct com- 
manded response values were essentially attained, with only slight changes 
noted for the control requirements. With regard to the 6-command maneuver, 
the coupled responses were negligible. Coupling was more pronounced for the 
other two commands but it was not appreciable. The Y“C°™"and maneuver exhib- 
ited a coupling of about 0.15 m/sec in speed. The largest amount of coupling 
was evident in the u-command maneuver where y approached 0.5°. This coupling 
would be much more pronounced for a large speed change, such as in high decel- 
ating approaches. As previously noted, however, the coupling effect could be 
eliminated by retrimming the flight-path angle. 

Effect of thrust .- The effect of a change in thrust coefficient on the 
decoupling process was not investigated. Data in reference 7 show that the 
effects of changes in thrust on the airplane stability and control derivatives 
(A and B matrices) are on the same order of magnitude as the effects of a 
changes; hence, similar effects on the decoupling process are indicated. 
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Summary of effects .- From the foregoing analysis of the effects of flight- 
condition changes, use of constant gains in a decoupled control system for cer- 
tain flight regimes (for example, during approach and landing) appears feasible. 
It has been shown that changes in the flight conditions do not cause excessive 
coupling in the response quantities (except for effects on Y during large 
speed changes) and have only minor effects on the control requirements. 


Effect of Airplane Control Lag 

The foregoing analysis assumed no control lag; that is, for the step com- 
mands, the control surfaces and engine thrust acted instantaneously. In any 
practical situation, control lag is always present in the flap and horizontal- 
tail servos as well as in engine thrust with respect to throttle movement. The 
effects of some typical values of lag are shown in the time histories presented 
in figures 14 to 17. The values used for and T3 were the same 

as those shown in table IV (a). 

Method of determining lag effects .- In order to determine the lag effects, 
lag was incorporated into the system as shown by the A and B matrices in 
table XI. For thrust lag effect (table XI (a)), fi^h considered as a state 

variable; that is, the actual (with lag) throttle or thrust control. Hence, 
the first column of the original B matrix is included in the fifth column of 
the A matrix. A first-order lag effect was assumed so that 


for time derivative (d/dt) , is the 

the thrust response time constant. The above 
in A and B in table XI. A value of 2.0 sec 
was used for ith* decoupling, the same G matrix was used as the one for 

no lag. The F matrix was the same but a fifth column of zeros was added. 

The procedure for investigating the effect of flap lag was the same as that 
shown for thrust lag. The third column of the original B matrix was included 
in the fifth column of the A matrix and first order lag was assumed with a time 
constant of 1.0 sec. 

For the horizontal tail, a second-order lag effect was assumed so that 


or 


'^th,! = f’^thD + D'Sth 


Kh 

6^h = - — + 
^th 


^th,l 

^th 


\rtiere D is the operator symbol 
input to the servo, and is 

equation appears as the last row 


■ o 

k6[,i = (p2 + D2CtO)t + a)t)6t 
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or 


v^ere K is merely a proportionality constant relating servo position ^t,l 
to control surface position 6^. The quantities 6^ ^t were con- 

sidered state variables and were incorporated into the A and B matrices as 
shown in table XI (b) . The second column of the original B matrix is included 
in the fifth column of the A matrix. As in the case of the other control 
lags, the original F and G matrices, calculated for ho lag, were used in 
determining the effect of horizontal-tail lag. The F matrix required addi- 
tional fifth and sixth columns of zeros. 


Effect on pitch-command maneuver .- The time histories in figure 14 show 
the, effect of lag in the different controls for a pitch-command maneuver. The 
horizontal-tail lag has a small time constant, is well damped, and as shown in 
figure 14(a), has no effect on the decoupling process. (This was also found 
to be true for the y- and u-command maneuvers.) The effects of lag in the other 
two controls are minor for a 0-command maneuver (figs. 14(b) and 14(c)). Except 
for the initial rates, the control requirements were not appreciably changed and 
the pitch-command value of 3° was essentially achieved. 


Effect of thrust lag .- The effect of thrust lag is shown in figures 15 
and 16, respectively, for a Y command and u conmand. Results are also 
included for and 0 to show the effect that these quantities 


'th 


have on decoupling with thrust lag present. The servo-input curves, as shown 
in figure 14, are not included in figures 15 and 16. The response quantities 
and control requirements for the no-lag condition are shown for comparison when 
lag is included. There are noticeable differences in the control requirements 
in all instances. The coupling effects do not appear to be excessive, except 
for the relatively large Au of close to -2 m/sec in figure 15(a). In several 
instances the lag has a large effect on the commanded response, with the most 
severe effect appearing in figure 16(a) where there is about a 70-percent over- 
shoot in Au. The y coupling in figure 16(b) would be much more pronounced 
in high decelerating approaches but could be corrected by re trimming the flight- 
path angle. 


Effect of flap lag .- The effect of lag in the flap for a Y~command maneuver 
and b-command maneuver is shown in figures 17 (a) and 17 (b) , respectively. Here 
again, the flap-servo input curves are not shown. The effects are similar to 
those shown in figures 15 and 16 for thrust lag, although not as pronounced. 

(As for the case of thrust lag (fig. 16(b)), large flap-lag effects on Y cou- 
pling would have to be trimmed out during high decelerating approaches.) Most 
notably, the changes in the control requirements are small, except for the 
initial rates. The flap rate in figure 17(a) of about 30°/sec during the 
Y-command maneuver (with lag) is well above the 5®/sec limit (table I) ; however, 
it applies to the extreme case of an instantaneous step-command input. Also, 
the correct command responses were attained in figures 17(a) and 17(b). It is 
of interest to note that the horizontal tail essentially takes over for the 
deficiency in the flap. 
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Suntmary of effects .- The effects of control lag on the decoupling process 
were more pronounced than the effects of other error sources investigated. The 
magnitude of the effects varied for the different controls and depended on the 
type of command maneuver. Pitch-command maneuvers were essentially unaffected 
by lag in any of the controls. Also, the effect of horizontal-tail lag was 
negligible for the different command maneuvers. Thrust-lag effect was the most 
critical with regard to coupling of responses and to attaining a desired command 
response. Except for the initial rates, the effects on the control requirements 
were not excessive. The coupling effects were generally moderate, with about 
0.5° in Y and 9, about 0.5 m/sec in u. Coupling in y would be much more 
pronounced for larger speed changes, but could be corrected by retrimming the 
flight-path angle. 


Inclusion of Lag in Decoupled Design 

The foregoing cinalysis clearly indicates that control lag could seriously 
affect decoupling. However, lag can be taken into account in the decoupling 
design, provided the lag characteristics of the controls are known. The pro- 
cedure used for including a first-order lag effect in the design of a decoupled 
control system is presented in appendix B. Example time histories with the con- 
trol lag included are shown in figures 18 to 20 for the three different command 
maneuvers. Perfect decoupling was achieved in all cases. In all cases the 
horizontal-tail servo time constant and the thrust time constant T^h were 

assumed to be 0.2 sec and 2.0 sec, respectively. Two different values were used 
for the flap servo time constant, 0.2 and 1.0 sec. Different values were used 
for the airplane dynamic response constants, as indicated in the figure titles. 
These airplane constants are noted only for the commanded response quantity 
inasmuch as the constants for the other response quantities were irrelevant* 
because these quantities were held to zero. 


Figure 18 shows exeimples of a Y-command maneuver for the two types of 
airplane dynamic response characteristics described in appendix B. (See 
eqs. (B14) and (B16).) As can be noted from appendix B, the results in fig- 
ure 18, as well as in figures 19 and 20, correspond to the case where Cgg ^ 

and ^ 0* The results in figure 18(a) are shown for two different val- 


ues of flap servo time constant Tf. The only difference in the time histories 
is in the flap servo input command. A good response in y is obtained, but it 
is obtained at the expense of an excessive initial flap rate of about 40°/sec. 
(See table I.) The high value of the throttle command is relatively insignifi- 
cant because the actual response value is very small. In figure* 18(b), a much 
lower initial flap rate (about 10°/sec) is required; however, the Y response 
is much slower because of the change in the airplane dynamic response character- 
istics. Even though the flap rate exceeds the limit of 5°/sec, it should be 
remembered that these time histories represent the extreme case of instantaneous 
step-command inputs. 


An example of a pitch-command maneuver is presented in figure 19 for the 
airplane pitch-response characteristics indicated in the title. The dynamic 
characteristics of the combined airplane- and control-response lags, as 
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described in appendix B, are included in the figure. The results show that an 
adequate pitch response is obtained. A high initial horizontal-tail rate of 
about 30°/sec is required; this is well within the limit value (table I) . 

Figure 20 presents various examples of a u-command maneuver. The different 
combinations used for the control- and speed-dynamic response characteristics 
are described in appendix B. (See eqs. (B14) and (B16) .) All the cases in 
figure 20 show that the actual response values of the flap and horizontal tail 
remain relatively close to the command values, whereas the opposite is true for 
the throttle owing to its higher t of 2.0 sec. Because the throttle is the 
primary control for speed in the decoupled design, this result indicates that 
the speed response ordinarily is affected appreciably if thrust lag is taken 
into account in the decoupled control system design. The of 0.2 sec used 

for the flap is unrealistically small; however, increasing it to larger values 
would affect only the flap-command value. 

Figure 20 (a) shows that a good speed response is obtained for the dynamic 
characteristics that were assumed. This response, however, was obtained at the 
expense of an excessive initial flap rate of 40°/sec. 

In comparison to figure 20(a), figures 20(b), (c) , and (d) show slower 
speed responses; however, the control rates are more reasonable.' In fig- 
ure 20(b), the speed response is not adequate; the actual throttle response 
value fails to reach the command value. In figure 20 (c) , the speed response 
time constant T 2 has been reduced to 0.5 sec. The speed response is better, 
but still too slow. The dynamic characteristics calculated for the combinations 
of Ti and are equivalent to the second-order system shown in the sub- 

title. Figure 20(d) is the same as figure 20(c), except that the damping ratio 
has been reduced by one-half. It can be seen that the speed response is greatly 
improved; however, about twice as much flap deflection is required for the 
maneuver. 

Summary of lag analysis results .- The foregoing results give only a brief 
look at the design possibilities which incorporate control lag in the decoupled 
system. Only one value for thrust lag and one value for horizontal-tail control 
lag were investigated together with several different airplane dynamic responses. 
There are countless ccanbinations of control an<3 airplane response characteristics 
that could be applied in the basic design of a decoupled system. The responses 
would necessarily have to be selected to match the control limits. 

The preliminary results presented in this paper indicate that, with proper 
selection of the airplane dynamic characteristics and appropriate adjustment 
of the control-lag characteristics, a decoupled control system can be designed 
to account for the adverse effects of control lag. For the 6- and Y-command 
cases, there appeared to be no problem in attaining adequate command response 
and reasonable control requirements. The u-command maneuver was the most crit- 
ical and required careful selection of the airplane dynamic characteristics. 
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CONCLUSIONS 


A Study has been conducted to determine the effect of various error sources 
on a decoupled control system. The analysis was performed in connection with 
the longitudinal-control system on a simulated externally blown jet-flap STOL 
aircraft. The system employed the throttle, horizontal tail, and flaps to decou- 
ple the forward velocity, pitch angle, and flight-path angle. 

It was shown that, despite the large errors employed in the study> the 
effects of the various errors on the decoupling process were generally minor in 
nature. Specific conclusions reached in regard to each of the error sources 
were: 


1. Effects of imperfect knowledge of the airplane characteristics (elements 
of the A and B matrices) were generally small. Decoupling was completely 
lost in only a few cases and this loss applied to a 50- to 100-percent error in 
the corresponding matrix element. As contrasted to B-matrix errors, the effects 
of A-matrix errors on decoupling as well as on the control requirements were 
essentially linearly related to the magnitude of the errors. This result indi- 
cates that designing adequate decoupled control systems requires better knowl- 
edge of the B matrix than of the A matrix. 

2. The dominant effect of state-variable measurement error occurred in the 
command response quantity in which the error existed; the control requirements 
and coupling in the other response quantities were not materially affected. 

The effect of measurement error was found to be nonlinear with respect to the 
magnitude of the error and (in the case of flight-path angle and speed) with 
respect to the nominal value of the pertinent airplane dynamic response charac- 
teristic defined in the A matrix. Because the error in the flight-path angle 
and speed-command response quantities depends on the magnitude of the nominal 
value of the A-matrix element, this error could be reduced if the A matrix 
were changed by altering the aircraft aerodynamic design. 

In general, the effect of measurement error on the command response would 
not be critical in a decoupled control system. The pilot would normally recog- 
nize this effect on the response he is commanding and could compensate for this 
error with an additional command input; only minor additional effects on the 
decoupling process would result. 

3. Except for flight-path angle coupling during large speed changes (for 
example, high decelerating approaches) , changes in the flight conditions did 
not cause excessive coupling in the response quantities and had only minor 
effects on the control requirements. Thus, use of constant gains in a decoupled 
control system for certain flight regimes appears feasible; consequently, the 
need for an onboard computer is avoided. The major effect of flight-condition 
changes occurred in the commanded response quantity; this effect, however, could 
be recognized by the pilot and could be easily controlled with an additional 
command input. 

4. The effects of control lag on the decoupling process were more pro- 
nounced than the effects of other error sources. Thrust-lag effect was the 
most critical with regard to coupled responses and with regard to attaining a 
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desired speed-command response, especially during high decelerating approaches. 
Control-lag effects on the control requirements were not excessive, except for 
the flap during flight-path-angle command maneuvers. 

5. The excessive flight-path-angle coupling caused by control lag and by 
changes in flight condition during speed-command maneuvers could be eliminated 
by retrimming the flight-path angle, either automatically or by the pilot. 

6. Provided the lag characteristics of the controls are known, lag could 

be incorporated in the design of the decoupled system. Such incorporation would 
permit complete decoupling; however, proper selection of the airplane dynamic 
characteristics and appropriate adjustments in the control-lag characteristics 
would be required. For pitch-command maneuvers there appeared to be no problem 
in attaining adequate command response and reasonable control requirements. 

The speed-command maneuver was the most critical and would normally require 
careful selection of the airplane dynamic characteristics. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 17, 1978 
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APPENDIX A 


DECOUPLED LONGITUDINAL CONTROLS 


The three longitudinal equations of motion were linearized as perturbations 
about an equilibrium condition in equations (1-59) of reference 13. These three 
equations can be solved simultaneously to give 


d2e' 


2 

UQ 


2pKyc2 


c I <39 ' 


UQ' 


dt 


Cma 


5 ^^ 

4y 


a' + C» + — — ^lu' 

. V “ J 




th 






'^«th ^ “IT ^ 


4y 


4y 


’f 4y 


da' UQ / c d0' , , 


(Al) 

(A2) 


du 

dt 


UQ 


2yc 


c^Q 




+ 1 Cx r + 

*^th 




4y 


^th + lCx5^ + 


%^Z6i 


4y 


K + |cx 5 + 


4y 


«f 


(A3) 


The terms C^ and Cg given in reference 13 were neglected. Also, sin 0 
^a q 

was assumed to equal 0 and cos 0 to equal 1 (0 is the angle between the hori- 
zon and X-equilibrium axis). The primed parameters in equations (Al) , (A2), 
and (A3) are perturbations from the equilibrium or trim conditions of the air- 
plane in nondimens ional form; that is. 


0' = 0 - 0Q 


(A4) 
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a' = a - Oq 


(A5) 


u - UQ 

u' = 

UQ 


(A6) 


and v^ere 


m 

P = (A7) 

PSc 



(A8) 


The weight and dimensional characteristics of the simulated airplane are 
presented in table I. The values used for the basic aerodynamic coefficients 
are given in table II. Constant coefficients were employed in the linearized 
longitudinal equations of motion corresponding to an angle of attack of 10°, 
a forward velocity of 30.48 m/sec, and a thrust coefficient of 1.87. 

The linearized longitudinal equations of motion can be written in state 
vector notation (after adding the fourth equation S = q) as 


X = Ax + Bu 


(A9) 


v^ere the state vector is 



(AlO) 


and the control vector is 


u = 




(All) 
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The airplane has a spoiler control, but it was not used in the study and, there- 
fore, was not included in equations (Al) to (A3). The control law used is given 
as 


u = Fx + Gv (A12) 

where v is the vector of commanded pilot inputs u^,, and that are 

to be controlled in a decoupled manner. The output equation 

y = Cx (A13) 

determines which state outputs are to be decoupled. 

Substitution of equation (A12) into equation (A9) yields the augmented 
equation 


X = (A + BF) X + BGv 

The Laplace transform of the result can be written as 


(A14) 


sx(s) - x(0) = (A + BF) x(s) + BGv(s) (A15) 


where the circumflex indicates the Laplace transform. Equation (A15) can be 
written as 




/N 


(sl - A - BF) x(s) = x(0) + BGv(s) 


(A16) 


Substituting the Laplace transform of equation (A13) into equation (A16) 
gives the transfer function between y and v as 

H = C(sl - A - BF)-1 bG (A17) 

For a completely decoupled system, the matrices F and G must be chosen 
so that the transfer function matrix H is diagonal and nonsingular. The com- 
puter program reported in reference 9 can be used to compute the desired F 
and G matrices. 
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Although the Gilbert program of reference 9 is based on a sophisticated 
algebraic theory of linear systems, the derivation of the decoupling control 
law (eq. (A12) ) can be simply explained by the method of Rekasius (ref. 1). 

The method is based on time differentiating each y-element until control quan- 
tities u appear explicitly in the expression for y; the associated v is 
then defined to give desirable transfer functions from v to y and to permit 
solving for the decoupling control law in the form of equation (A12) . For pur- 
poses of clarity, the primes are dropped from the response and control variables. 

For speed control only, let 


Yl = u ; Yl = u 


(A18) 


From equation (A3) and table III , then 


yi = 94,1® + 94,3“ + 94^4U + b4^2<^th + b4^2<^t + ^ 4 ^ 36 f 


(A19) 


(Note that eq. (A19) contains required control quantities.) Adding first-order 
dynamics yields 


yi + — = 34^20 + 34^30 + (a 4^4 + — ju + b4^i6th + b4^2^t *^4,3<^f = (A20) 


Similarly, for pitch control only, let 


Y2 = 6 ; y2 = ® ; Y2 = ® 


(A21) 


From equation (Al) and table III, then 


Y2 = 92, 2q + 92,3“ + 92,4“ + b2,i6th + b2,2<St + b2,3fif 


(A22) 


(Eq. (A22) contains control quantities.) Adding second-order dynamics yields 


Y2 + 2CUnq + 0)^0 = 0)^0 + 


( 92,2 + 2Ca>n)q + 32,3a + 32,43 + b2,i6th 


+ b2,2^t ‘^2,3'^f * ^2 


(A23) 
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Finally, for flight-path-angle control only, let 

y3=0-a;y3=e-a (A24) 

From equation (A2) and table III (and noting that 6 = q) , then 

y3 = (1 - a3,2)q “ 33,3°^ “ 33,4^ - b3,i6th - b3,2<5t - b 3 , 36 f (A25) 

(Bq. (A25) contains control quantities.) Adding first-order dynamics yields 


a 

y3 + — = (1 - 33 2)q 
•^3 




a - 33,4“ - b3. 


I'^th 


- b 3 , 2 fit - b3,36f = V3 


(A26) 


Equations (A20), (A23), and (A26) can be collected to form 


l\ 

VI = 34,10 + Oq + a4,3« + [a4,4 ~j'^ * *^ 4 ,l^th + b 4 , 26 t + b4,36f 


V2 - “n® + (32,2 + 2C(Jn)q + 32,3“ + 32, 4 U + b2,i6th + ^2,2^t + b2,3<5f / (A27) 


V3 - 00 + (1 - a3,2)q - (33,3 - — )a - 33, 4U - b3,i6th - b 3 , 26 t " b3,3<Sf 

J 


Equations (A27) can be written in matrix form as 


V = F'x + G'u 


(A28) 


or 


G = -[g']"1 F'x + [G']-1 


(A29) 
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Providing G' is not singular, then 


u 4 Fx + Gv 


(A30) 


where F = -[G']“^F' and G = [G']“^ are the feedback and feedforward gain 
matrices required for decoupling. The values of T3, X,, and are 

free to be selected for any desired response. 
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DECOUPLED LONGITUDINAL CONTROLS WITH LAG 

In this section the longitudinal equations of motion given in appendix A 
are altered in order to take control lag_into account in the decoupling design 
process by defining the control vector u as 


u 


^t,l 


(Bl) 


v^ere the subscript 1 refers to servo control. In the design procedure, lags 
are considered first order by letting 


(TthD + l)<Sth = ^th,l 


or 



•^th 


(TfcD + l)6t = <5t,l 


or 




'■ . 5f «f,l 

(TfD + l)6f = <5f 1 or 5f + 

' Tf Tf 


(B2) 


(B3) 


(B4) 


where D is the operator symbol for the derivative with respect to time and 
6th r *^t» ^f become state variables. 

Equations (B2) to (B4) , combined with the longitudinal equations of motion 
of appendix A can be written in matrix form as 

« 

X = Ax + Bu (B5) 


27 



where 


28 


B 


(B6) 


in table III (a) is 


u 

<Sth 



0 

0 

0 

0 

0.2250 

-0.1460 

-2.380 

0.1487 

-0.6400 

-0.1850 

-0.0676 

-0.1712 

-0.1018 

0.1047 

-0.01406 

-0.1190 

0 

-1/Tth 

0 

0 

0 

0 

-lAt 

0 

0 

0 

0 

-1/Tf 


(B7) 
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and matrix B is 


0 

q 

a 

B = u 


5 th 


°th,l 

0 

0 

0 

0 

lAth 

0 

0 


°t,l 

0 

0 

0 

0 

0 

lAt 

0 


0 

0 

0 

0 

0 

0 

lAf 


(B8) 


^In the A matrix, the C2g and ^ 0 case is used.j 


The procedure for determining the F and G matrices required for decou- 
pling is similar to that shovm in appendix A. There is a difference in the 
required number of differentiations of y, as shown below. Here again, the 
primes are dropped from the response and control variables for purposes of. 
clarity. 

For speed control only, let 


Y1 = u ; yi = u 


(B9) 


thus, from equations (B7) and (B8) 


yi = 34, l9 + 34,3“ + 34, 4U + 34,56th + a4,e5t + 34,76f (BIO) 

It is evident, however, that equation (BIO) does not contain a servo-control 
variable which is required. Therefore yx must be differentiated again 

yi = ii (Bll) 

so that 

yi = 34, x6 + 34,3a + a4,4U + 34,56th + 34,6^t + 34,75^ (B12) 
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Substituting for the differentials 


yi = + 34 , 3 ( 33 , 2 '! + 33^3tt + a3^4U + a3^56th + 33,6*^t + 33^76f) 


+ a4^4(a4^10 + a4^3tt + a4^4U + a4^56th + 34^65t + 34^76f) 


* 34,7(- ^ ^ «f.l) 


(B13) 


(Eq. (B13) contains a servo-control variable which is required.) 

As for including dynamics in equation (B13) , the results indicate second- 
order dynamics for the u response; therefore, the u response is modeled as 


1 1 \ ' 1 . 1 \ 1 

Is + — s + or s^ + s — + + 

Tl/V Tth/ Wl Tth/ TiTth 


(B14) 


The speed response would then have the engine response lag superimposed oh the 
airplane response lag. Adding the dynamic response to equation (B13) yields 


'1 1 \ (I 

ii + ul — + 1 + u| 1 = vi 


V^l Tth; 


V'^l'^thy 


= 0 


1 1 

a4,4 + |— + — 
^1 ^th, 


34,1 + <!(34,1 + 34,3) 


+ a 


1 1 

33,3 + 34,4 + I— + — 
^1 ^th, 


34,3 


+ u 


, 1/11 

(34,4) + 34,3S3,4 + — + a4,4| — + 


T^l^th 


'^th, 


(Equation (B15) continued on next page) 
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+ fith\a4, 3^3,5 + ^4,5 


1/1 1 
34 4 + — + 

'^th \Tl Tth 




+ 'St<(a4,3a3,6 + 


1 1 1 

34,4 + — 

Tt \ti Tth 



+ 6 


f'^4,333,7 + 34, 7| 


1 /I . 1 

34,4 + — + — 

Tf \'^1 Tth 



^ ^ * «t,l(a4,6 * «£.l(=4,7 


(B15) 

Instead of using the expression given by equation (B14) , the second-order 
dynamics for the u response could be represented by 

2 

u + 2^]^o)]^u + (O^u (B16) 


1 1 

and the terms 1 — + — - 

Tth; 


and 


TlTth 


in equation (B15) would be replaced by 


( 2 ^l 0 )l) and 0 )i, respectively. 

Similarly, for pitch control only, let 


■y '2 = q (B17) 

so that, from equations (B7) and (B8) 

• • • 

‘y '2 = 32,2^ + 32,3& + 32, 4U + 32 , 56 th + 32,e6t + 32,76f 
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Substituting for the differentials, 

y '2 = 32, 2 (32,25 + 32,3® + 32, 4U + a2,5<Sth + 32,6'^t + 32,7^f) 

' + ^ 2 , 3 (33 , 25 + 33,3a + 33, 4U + 33,56th + 33 , 6 <St + 33 , 76 f) 

" + 32,4(34,16 + 34,3a + 34,'4U + 34,56th + 34,e6t + 34,765) 


.+ 





+ 




(B18) 


(Eq. (B18) contains servo-control variables.) 


As for including dynamics in equation (B18) , the results indicate third- 
order dynamics for the pitch response; therefore, it appears logical to assume 
for the response 



+ 2C0)„s 




(B19) 


The lag in the horizontal tail would then be superimposed on the airplane 
response lag. Adding the dynamic response to equation (B18) yields 


a + o I 2C0X, 


L 


\ 


1 \ / 2 1 \ A 2 \ 

+ — +0 0 )., + — 2C0)., 1 + 91 — uo.. 1 = VO 

V V' 7 7 


(B20) 


or 


q + qP + qQ + 0 R = V2 = 9 ( 32 , 434,1 + R) + 5 ( 32 , 2 )^ + ^ 2 , 2 ^ 2, 2 + Q + 32 , 2? 


“[32,3(32,2 + 33,3 + P) +.32,434,3^ 


(Equation (B21) continued on next page) 
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APPENDIX B 


“1^32,4(32, 


2 + 34,4 + P) + 32,333,4 


+ 6 


th 


+ 6i 



J 

- 

f 1 \ 

32 , 333,5 + 32,434,5 + 32,51 

32,2 - ^ + P| 


^ ^th /. 

n 


32,333,6 + 32,434,6 + 32,6 32,2 + P 

'^t 


+ 6^ 


32,333,7 + 32,434,7 + 32,7(32,2 “ ~ 


+ ^th,l(32,5 — j + «t,lf32,6 + <Sf,ifa2,7 — 


(B21) 


where 


P = 2CWn + — 


2 1 

Q = Wjj + — 2^c0n 
■^t 


T, ^ 2 

R = — Wn 


The equation for V3, representing y control only, can be obtained by 
letting y3 = 0 - a and following the procedure for V3. (See eqs. (B9) 
to (B15).) The equations for vi, V2» and V3 can be combined in like man- 
ner to, equation (A27) and solved for the gain matrices P and G as in 
equation (A30). 
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TABLE I.- WEIGHT AND DIMENSIONAL CHARACTERISTICS OF SIMULATED AIRCRAFT 


Weight, N . . . 

Wing area, 

Wing span, m . 

Mean aerodynamic chord, m 


Center-of-gravity location, 
ly, kg-m2 ......... 


percent 


c 


Maximum control-surface deflections; 

<3eg 

6f, deg 


Maximum control- surface deflection rates 

6^, deg/sec . . 

5f , deg/sec 


245 096 
78 
24 

. 3.58 

. 40 

334 637 


±10 
0 to 90 


50 

5 
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TABLE II.- BASIC AERODYNAMIC CHARACTERISTICS USED IN STUDY 


[values in parentheses used for extra cases] 


Parameter 

Value for - 

a = 10° 

a = 5° 

Cvj(*C2) 

-5.71 

-5.10 

% 

-6.567 

-7.20 

% 

2. 81 

1. 97 

Sx 

-1. 424 

-1. 80 


0 

. \ 

0 


-13.14 

-13.70 

Cx 

, 

0 

0 

C 7 

-11. 42 

-10. 20 

Cx 

*u 

-1. 82 

-1. 40 


0 

0 


-32.86 

' -34.30 

C7 

<5t 

-1. 209 

-1.157 


-. 252 

-. 0688 

^6t 

-5.272 

-5.56 





w, \-3. j; 

U 

^5th 

0, (-.468) 

0 

"*5th , 

1. 87 

1.87 


-3.06 

-3.40 


-2.126 

-2^225 

^6f 

.195 

.126 
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TABLE in.- VALUES OF COEFFICIENTS FOR A AND B MATRICES USED IN STUDY 
[values in parentheses used for extra cases] 





(a) 

a = 10° 









Throttle 

Horizontal- 

Flap 


e 

q a , 

u 



tail 


• 

9 

i 

0 

1.0 0 

0 

• 

0 

0 

0 

0 

• 1 

q 

0 

-1.23 -0.52 

0.225 

• 

q 

0 



A = . 





(-0.146) 

-2.38 

0.1487 

a 

0 

1.0 -0.368 

-o:64 

B = 




• 




• 

a 

0 



u 

-0.3195 

0 0.157 

-0.1018 


(-0.185) 

-0.0676 

-0.1712 





• 

u 

0.1047 

-0.01406 

-0.1190 


(b) a = 5° 


Throttle Horizontal- Flap 

0 q ct u tail 


• 

" 



- 

• 

“ 


- 

0 

0 

1.0 . 

0 

0 

0 

0 

0 

0 

• 

q 

A = . 

0 

-1.285 

-0.6725 

0.209 

• 

q 

B = . 

0 

-2.515 

0.127 

a 

0 

1.0 

-0.403 

-0.57 

CX 

0 

-0.0646 

-0.19 

• 

u 

-0.286 

0 

0.1105 

-0.0785 

• 

u 

0.1047 

-0.00383 

-0.1247 
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TABLE IV.- TYPICAL GAINS TOR DECOUPLING LONGITUDINAL 


AIRPLANE MOTIONS U, 0, AND y 


(a) Standard case 


= 2 

rad/sec 

= 1.0 

sec 

5 = 0.7 


T3 = 1.0 

sec 

-3.99115 

-0.20247 

2.69253 

-12.78503 

1.28404 

0.64378 

0.01187 

-0.13568 

-6.34813 

-0.25420 

3.68691 

-3.68474 


9.55110 

0.12896 

6.52689 

0 

-0.41005 

0.35616 

0 

0.16191 

5.70049 


(b) Ctt and ^ 0 case 

^th ^th 


fiJn = 2 rad/sec = 1.0 sec 

C = 0-7 T3 = 1,0 sec 


F 


-1. 82356 

1.51337 

-4.46814 


-0.09250 

0.65541 

-0.15885 


1. 23023 
-0.14285 
2. 41860 


-5.84152 
0. 59894 
2.33758 


G 


4.36393 

-0.54880 

-4.49900 


0.05892 

-0.41746 

0.10117 


2.98216 

-0.01887 

2.62603 



TABLE IV.- Concluded 


(c) Case with lag included in the decoupling design 

= 2 rad/sec = 0. 2 sec 

C = 0. 7 Tf = 1. 0 sec 

= 0. 5 sec = 2. 0 sec 

T3 = 1. 0 sec 


0.22497 

-0.23249 

0.76821 

1.58269 

-0.30732 

-0.04432 

0.05331 

-11.77187 

-0.62304 

-1. 88429 

-0.40576 

1.62150 

-3. 03945 

1. 39079 

0.19593 

2.54177 

-2.64842 

-7.55244 

1.11315 

0.11508 

-1.04370 



-0.10976 

-0.08349 

-0.00377 




G = 

8.72786 

0.11785 

5.96432 





-4.49900 

0.10117 

2.62603 



















TABLE VI.- EFFECT OF UNCERTAINTY IN B-MATRIX COEFFICIENTS ON DECOUPLING 


6 command 


q 0 



“ ^X6 


‘Quantities have minor effect on 
element. 


y command 




*9 undershoot, 
12 percent. 


y undershoot, 
28 percent. 


u command 


■ 

■ 

■ 

■ 

■ 

■ 


H 

a 

■ 


*20 percent u 
overshoot for 
20 percent u under 
shoot for b 4 ^ 3 . 




I I Value cannot change 

No significant effect on decoupling 
00B Decoupling lost for 9, y, and u, respectively 
000 Partial decoupling lost for 9, y, and u, respectively 
n The terms "overshoot" and "undershoot" refer to steady-state error. 






TABLE VII.- SUM^Y OF EFFECTS OF INDIVIDUAL MEASUREMENT ERRORS 

ON DECOUPLING PROCESS 

[Blank blocks indicate no significant effect] 


Measurement 
error in - 

Response 

commcinded 

Effect on - 

0 

Y 

u 

0 

0 

Large 

Small 

Small 

Y 




u 




• 

0 

0 




Y 




u 




a 

0 

Negligible 

■ 

Small 

Small 

Y 

Small 

Large 

Small 

u 




u 

0 




Y 




u 



Negligible 

Small 

Large 
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TABLE VIII.- COMPARISON OP EFFECTS OP COMBINED MEASURH4ENT ERRORS 


ON DECOUPLING PROCESS 

[All control differences from nominal are small cuid less than those shown 
in time history in figure 8 for 40-percent error in 9] 


(a) 3° command in 0 


Condition 

Maximum response in - 

9, deg 

y, deg 

u, m/sec 

• 

Positive 20-percent error in 0, d, a, and u 

2.5 

-0.15 

0.035 

• 

Positive 20-percent error in 6, 9, and u; 

negative 20-percent error in a 

2.4 

-.90 

.110 

Positive 20-percent error in 0 only 

2.5 

-.45 

.075 

Negative 20-percent error in 0 only 

, — ^ 1 

3.7 

^ 

.70 

-.120 


(b) 6° command in y 


Condition 

Maximum response in - 

y, deg 

9, deg 

u, m/sec 

• 

Positive 20-percent error in 0, 9, a, and u 

5.3 

-0.11 

0.070 

• 

Positive 20-percent error in 0, 9, and u; 

6.7 

.14 

-.080 

negative 20-percent error in a 




Positive 20-percent error in a only 

5.3 

-.13 

.085 

Negative 20-percent error in a only 

6.7 

.16 

-.100 
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TABLE IX.- EFFECT OF SPEED ON A AND B MATRICES 

“ [Multiplication factors for effects in A and 
B matrices due to speed changes. (Blank 
blocks denote no change)] 



A 

B 

0' 

q' 

a' 

u' 

&th 

6t 

6f 

d9 '/dt 







■■ 

d20 Vdt2 


u/UQ 

(u/uo)2 

U/UQ 

(u/uq)2 

(u/uq) 2 

(u/uq)^ 

da '/dt 



u/uq 


U/UQ 

U/Uq 

u/uq 

du ' /dt 

U/UQ 


U/UQ 

■ 

U/UQ 



u/Uq 

u/uq 


44 


























TABLE X.- EFFECT OF CHANGE IN FLIGHT CONDITION ON DECOUPLING PROCESS 


[Blank blocks denote no effect] 



Speed increased 
by 50 percent® 

Speed increased 
by 100 percent 

Speed decreased 
by 50 percent 

Maximum response due to 1.5 nv^sec command in u 

u, m/sec 

1. 54 (2.44) 

• 

1. 58 

1. 35 (sluggish) 

0 , deg 

-0.055 (0.017) 

0 

• 

0 

1 

0.15 

Y» «3eg 

-0.63 (-0.20) 

- 1.0 

1.4 

Maximum response due to 3° command iri 6 

u, m/sec 


0 . 05 


9, deg 

3.0 (3.0) 

3.0 

3.5 

Y, deg 


-0.025 


Maximum response due to 6 ° command in Y 

u, m/sec 




0 , deg 




Y» deg 

6.0 ( 6 . 0 ) 

6.0 

6. 0 (sluggish) 


^Values in parentheses from ^ 25 ^^ ^5 h ^ ^ case. 
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TABLE XI.- A AND 


B 


MATRICES WITH EFFECT OF LAG INCLUDED IN SYSTEM 


(a) Lag in thrust® 



e 

q 

a 

u 

«th 

• 

e 

0 

1.0 

0 

0 

0 

• 





0 

q 

0 

-1.23 

-0.52 

0.225 

(-0.146) 

• 





0 

> 

II 

R 

0 

1.0 

-0.368 

1 

o 

• 

(-0.185) 

• 

U 

-0.3195 

0 

0.157 

-0.1018 

0.1047 

Kh 

0 

0 

0 

0 

-0.5 



*^th,l 

<St 




• 

e 

0 

0 

0 



• 

q 

0 

00 

« 

CM 

1 

0.1487 



• 

B = a 

0 

-0.0676 

-0.1712 



• 

u 

0 

-0.01406 

-0.1190 



«th 

0.5 

0 

0 



^Values in parentheses from C^^j ^ and C„,g ^ 0 


case 



TABLE XI.- Concluded 


(b) Lag in horizontal- tail servo 


-1. 23 -0. 52 


0. 225 


-2.38 


-0.368 


-0.64 


-0.0676 


0.157 


-0.1018 -0.01406 


*^th *^t,l 


-1046 -95.4 


q 0 


0.1487 


a 0 


0 -0.1712 


u 0.1047 0 -0.119 


6 ^. 0 
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Figure 2.- Comparison of pitch maneuvers performed by unaugmented and 
decoupled control systems. (Nominal a = 10°. ) 





2 3 

Time, sec 

3° pitch command 


4 


5 


*f errors in A-matrix coefficients on 
Lons and commanded responses. 




















Nominal (controlling 7 or a) 

1+0 ^ error In 9 (controlling 7 ) 
1*0 ^ error in 9 (controlling a) 









Measurement error In command 


0 



Au, 

m/sec 


.2 


Ar, 

deg 


- 21 — 



Figure 11.- Effect of 






(b) 6° flight-path-angle command 


Figure 13.- Continued. 

















Au, 

m/ sec 



A/, 

deg 


A9, 

deg 


^thJ 

percent 


A5f, 

deg 


ABt, 

deg 


0 



I I I I I I 

0 1 2 5 4 5 


Time, sec 

(c) Lag in engine response = 2.0 sec). 

Figure 14.- Concluded. 
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Figure 15.- Effect of engine-response lag (T^j^ = 2.0 sec) for 
flight-path-angle command of 6°. 









With lag 


All, 

m/ sec 



A9, 

^eg 



-5 




(b) Cgc and Cn« 0. 
^th ^th 


Figure 16.- Concluded. 
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0 1 2 5 4 

Time, sec 


(a) Flight-path-angle command of 6°. 

Figure 17.- Effect of flap-servo lag (tf =1.0 sec) 




Input to servo 1.0 sec) 


— Input to servo (T^ = 0.2 sec) 


Actual response (with lag) 














(a) 

Figure 20.- Speed-command ma 
into design of decoupled 
T^h = 2.0 sec). 
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